Plasma Suppression of Large Scale Structure Formation in the Universe 
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We point out that during the reionization epoch of the cosmic history, the plasma collective effect 
among the ordinary matter would suppress the large scale structure formation. The imperfect Debye 
shielding at finite temperature would induce a residual long-range electrostatic potential which, 
working together with the baryon thermal pressure, would counter the gravitational collapse. As 
a result the effective Jean's length, A./, is increased by a factor, Xj/Xj = a/8/5, relative to the 
conventional one. For scales smaller than the effective Jean's scale the plasma would oscillate at the 
ion-acoustic frequency. The modes that would be influenced by this effect depend on the starting 
time and the initial temperature of reionization, but roughly lie in the range 0.5/iMpc -1 < k, which 
corresponds to the region of the Lyman-a forest from the inter-galactic medium. We predict that in 
the linear regime of density-contrast growth, the plasma suppression of the matter power spectrum 
would approach f — (fldm/^lm)' 2 ~ 1 — (5/6) 2 ~ 30%. 
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The matter power spectrum of the large scale struc- 
ture, alongside with the cosmic microwave background 
(CMB) fluctuations, provides important information on 
the composition and the evolution history of the uni- 
verse. Unlike the CMB fluctuations, which was frozen 
since the recombination time, the matter power spec- 
trum continued to evolve until now, during which the 
universe has been essentially in a reionization, or plasma, 
state. Plasmas are known to support numerous electron 
and ion collective oscillation modes. Although in the 
pre-recombination era the universe was also in a plasma 
state, the extreme dominance of radiation pressure over 
that of matter (by a factor 10 10 ) before decoupling ren- 
ders the contribution of (ep) plasma oscillations totally 
negligible, fn the post-decoupling era, the CMB pressure 
became negligible. Thus during the reionization epochil] 
one might expect collective plasma effects to exhibit in- 
fluence on the structure formation. 

A perfectly uniform plasma is charge-neutral and has 
no long-range electrostatic potential. The inevitable 
thermal fluctuations of its electron density at finite tem- 
perature would, however, induce a non- vanishing resid- 
ual electrostatic potential among the ions that defies the 
Debye shielding. Such a residual electric force that ex- 
erts an additional pressure on the ions is non-collisional 
in nature. When combined with the collision-induced 
thermal pressure among them, these pressures can sup- 
port ion-acoustic oscillations, a phenomenon well-known 
in plasma physics. [2|, y] 

In this Letter we point out that the inclusion of such 
plasma collective effects would result in an increase of 
the sound speed and therefore an increase of the Jean's 
scale by a factor Xj/Xj = a/8/5. The modes that would 



be influenced by this effect depend on the starting time 
of reioniation and its initial temperature, but should 
roughly lie in the range 0.5/iMpc -1 < k, which corre- 
sponds to the region of the Lyman-a forest. This plasma 
suppression mechanism is effective to both linear and 
nonlinear regimes of density perturbations. We predict 
that, the maximum suppression of the power spectrum in 
the linear regime should approach ~ 1 — (f2d m /f2 m ) 2 ~ 
1 - (5/6) 2 ~ 30%. 

When the electrostatic potential is included, the stan- 
dard Einstein-Boltzmann equation in the Newtonian 
limit is now extended into a Maxwell-Einstein-Boltzmann 
(MEB) equation. Let us designate the baryon perturba- 
tion as A e Spb/pb. In the electrostatic and the New- 
tonian limits and with the density perturbation A C 1, 
the MEB equation reads, 



A + 2HA 
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where the electrostatic potential satisfies the Poisson 
equation: V 2 e m = 47rep e . The source of gravity con- 
tains both dark matter and baryons, i.e., p = pd m + Pb- 
We assume that the reionization epoch began at a cer- 
tain time a*. The adiabatic condition dictates that 
the initial perturbations are identical between the dark 
and the baryonic matters, i.e., 5pdm( a *)/pdm( a *) = 
Sp b (a*)/pb(a*) = 5p(a*)/p(a*) = A*. After that, dark 
matter and baryonic matter would still evolve hand-in- 
hand and grow linearly for scales larger than the Jean's 
scale, but would evolve separately for scales smaller than 
that, ft is this latter situation which we are focusing for 
our plasma effect. Thermodynamics governs that 
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On the other hand, the electrostatic potential at a given 
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temperature is given by 
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The approximation results from the fact that the typi- 
cal cosmic scale of interest is much lar ger than the De- 
bye length, i.e., A ~ 1/k > X D = yJ^ e k B T e /ATTe 2 p e . 
We have further assumed that 8p e /p e ~ Spb/pb = A. 
This is because globally the fast-moving electrons tend 
to rapidly readjust themselves to follow the baryon den- 
sity perturbations. At first sight, it may seem counter- 
intuitive that the charge neutrality condition would al- 
low for a non-vanishing electrostatic potential. In re- 
ality, due to the inevitable thermal fluctuations of the 
much lighter electrons the plasma is never perfectly neu- 
tral. Balance of forces on electrons requires that the 
electron number density satisfies the Boltzmann's re- 
lation: n e = p e /m e = n e0 exp(-e4> em /k B T e ). Thus 
at finite temperature the Debye shielding is never per- 
fect, and there is always a residual electrostatic poten- 
tial at the level of e<p ern ~ kBT e 5p e / p e in the lead- 
ing order of Taylor expansion. This, that one assumes 
Spe/Pe = Spb/pb and V • E = ~S/ 2 4> em ^ at the same 
time, is the so-called quasineutrality, or plasma approxi- 
mation in plasma physics 0. 

Fourier transforming Eq.JT]) leads to 



where 
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is the 'sound speed' of the ion-acoustic waves. We see 
that the main change from the ordinary acoustic wave to 
the ion acoustic wave is that there is an extra contribu- 
tion to the sound speed from the electron temperature. 
The fast- moving electrons are isothermal. So j e = 1. 
But for protons, jb = 5/3. In the case where electrons 
and ions are in thermal equilibrium, i.e., T e = Tb, the net 
change is an increase o f th e sound speed from the ordi- 
nary acoustic wave, to that of the ion-acoustic 
wave, v/1 + 5/3. Both waves, however, are constant- 
velocity pressure waves that satisfy linear dispersion re- 
lations in the long wavelength limit. 

So far the temperature and the on-set of the reioniza- 
tion are still unresolved. In the literature (e.g. 0,IEE|)> 
the reionization temperature ranges from 10 3 K to 10 7 K. 
As for the on-set of reionization, the observation of the 
Gunn-Peterson trough in the spectra of high redshift 
quasars indicates that the reionization process was com- 
pleted by z* « 6 @, 11, while WMAP determines that 
z* = 17 ± 5 0, For simplicity and without compro- 
mising qualitative features of our effect, we model the 
reionization as an abrupt transition at a given instant 
a* = 1/(1 + z*) where the plasma temperature raised 
instantly to T b * homogeneously throughout the universe. 



To solve Eq.® for a given fc-mode, it is essential to rec- 
ognize its relationship with the initial and the final Jean's 
scales at a* and a, respectively. Since the universe ex- 
pands adiabatically during the matter-dominant epoch, 
the large-scale baryon/plasma temperature satisfies the 
equation of state, Xf,V 7i>_1 = const., where jb = 5/3 and 
V cx 1/p oc a 3 . Therefore Tb oc c 2 oc a~ 2 . Thus the 
post-reionization (a > a*) Jean's scale evolves as 



kj(a) = 



^4:TTGp(a) _ /3vr m b Gp*\ 1 / 2 



,(a)/a 



T* 



'a*a > k*j . (6) 



There are three regimes of /c-modes that concern us. For 
k < kj, the thermal and the plasma pressures are negli- 
gible from the outset. So the baryon perturbation grows 
linearly in the same way as that of the dark matter: 
6b oc 8dm oc a. The decay mode, Sudm) a~ 3 / 2 , fades 
away in time and can be ignored. 

For k > k*j, care must be taken regarding fc's further 
relationship with the Jean's scale at the time of interest. 
If k > kj(a) > k*j, then the combined pressure dominates 
over the gravity throughout the period (a* , a) . Thus the 
right-hand-side of Eq. ^ is negligible. Changing variable 
from t to a, with a oc t 2 / 3, in the matter-dominant era, 
we find 
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Let us introduce a parameter uj r which satisfies the rela- 
tion, k 2 c 2 /H 2 = u 2 a~ 2 , such that ui r is independent of 
a. Since kj oc a 1 / 2 , we can reexpress it as 
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It can be shown, through one more change of variable 
to x = a~ 2 , that the solution to Eq.® is oscillatory: 
Afc oc exp(i2cj r / v / a). Physically uj r /y/a represents the 
phase of the ripples of the density contrast in fc-space due 
to the ion-acoustic oscillation. Such an oscillation begins 
at the time of reionization, a*. Although in reality c s k 
is the physical ion-acoustic oscillation frequency, uj r can 
be viewed as an effective ion-acoustic frequency since it 
relates to c s k straight-forwardly. 

Ion-acoustic oscillations are known to suffer Landau 
damping 2], which was not included in Eq.® a pri- 
ori. To include this effect, we invoke a simple empirical 
formula 2| for the Landau damping dispersion relation 
that relates the real and the imaginary parts of the ion- 
acoustic oscillation frequency, 



^~i.i(Ti) v \-(Te/T b y 

LU r V Tb J 



1 < T e /T b < 10 . (9) 



Since the electron and proton temperatures in the reion- 
ization epoch are roughly equal, we find that LUi/uj r ~ 0.4 
in our case. This means that within ~ 0.4 of an oscil- 
lation, the ion-acoustic wave would be Landau-damped 
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to one e-folding of its initial amplitude. Now we replace 
the cu r in the exponent of the solution to Eq.([7]) by the 
complex frequency which satisfies the above dispersion 
relation. Imposing the initial condition at a*, we arrive 
at the perturbation for the fc-mode whose scale has been 
below the Jean's scale all the time from a* to a, i.e., 
k > kj(a) > k*j, 



A fc ~ ^{-o-s+a^v . 



(10) 



where r\ = r\(a* , a) = Wlja* — Wlja. 

For the intermediate regime fc} < k < kj , the fc-mode 
must have exited the Jean's scale at an intermediate time 
a c : k — kj(a c ), where a* < a c < a. Such a mode would 
oscillate for a period Aa = a c — a*. Then it would stop 
the oscillation at a = a c and resume its gravitational 
collapse during the subsequent interval a — a c . Neverthe- 
less, such a resumed growth of perturbation is delayed, 
because of the ion-acoustic oscillations during the earlier 
time interval, by the same amount 



Aa = a r — a* = 



("J 



Thus the evolution of such a mode kj < k < kj(a) should 
become 
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Here we have ignored the fact that upon crossing the 
Jean's scale the fc-mode may still carry a non-vanishing 
oscillation amplitude. We expect that such remaining 
amplitude should be small due to Landau damping. 

Our aim is to identify the imprints of this plasma ef- 
fect in the large scale structure formation through obser- 
vations. So we look for its modification to the matter 
power spectrum [Til]. 



P(k,a) = 2n 2 6 2 H 



Tjn+3 
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(13) 



Here 6jj = 1.9 x 10~ 5 is the amplitude of primordial per- 
turbations, Ho — 100ft.kms _1 Mpc _1 the Hubble param- 
eter at present, T(k) the transfer function, and D{a) the 
growth function. We limit our investigation in the lin- 
ear regime of power spectrum growth. Due to the onset 
of the plasma effect, the conventional growth function, 
D(a), in Eq. lfTB"]) has to be modified. Let us denote the 
modified growth function as D(k, a). Clearly, the plasma 
suppression would only affect the regime where k > kj, 
which is what we will concentrate below. To track the 
plasma effect on the baryons for a > a*, we write 



D[k,a) = D dm (a)^ + D b (a) 



£M a ) — ~FT~ + D b (a )-T77^-. (14) 
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Invoking the initial condition, Dd m {a*) = Db{a*) — 
D(a*), and under the conventional notion of linear 
growth, D(a) — D(a*)(a/a*), we arrive at the relative 
change of the growth function as 



D{k,a) 



a* A k 
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Thus for a given k-mode with k > kj, the plasma- 
suppressed matter power spectrum relative to that of the 
conventional approach reads, for k > kj(a) > kj, 
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On the other hand we have, for fc} < k < kj{a) 
D(k,a) 2 _ {£ldm\ 2 U , n ( ^6 a 
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FIG. 1: Suppression of the matter power spectrum, a* — 0.1 
and T 6 * = 5 x 10 6 K. The red, purple and blue curves represent 
the suppressions at a = 1/6,1/3 and 1, respectively. The 
three regimes of k can be clearly recognized. 

To appreciate what this implies in cosmology, we look 
for the physical scales relevant to the plasma suppres- 
sion effect. Let us assume that the reionization occurred 
at 1 + z* = 10, or a* = 0.1, with the initial plasma 
temperature T* ~ T b * ~ 5 x 10 6 K. This corresponds 
to an ion-acoustic wave velocity of c* ~ 7 x 10 _4 c and 
kj ~ 0.63/iMpc _1 , and therefore the Jean's length, A} = 
2w/kj ~ 1.0/i _1 Mpc. The Jean's scale at present would 
be kj(a = 1) - 2.0/iMpc~ 1 , or Xj(a = 1) ~ S.l/i^Mpc. 
Figure [T] plots the suppression factor, \D(k, a)/D(a)\ 2 , at 
different times: a =1/6, 1/3 and 1, based on the ACDM 
model with = 0.7, O m = 0.3, = 0.25, Oj, = 
0.05, a* = 0.1, a = 1 and T* = 5 x 10 6 K. The kinks of 
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each curve are associated with the transitions at kj and 
kj(a c ), respectively. These kinks are the artifact of our 
approximation, which should be smeared in reality. The 
three regimes of k, namely, k < k*j, kj(a = 1) > k > kj, 
and k > kj(a = 1) > kj, can be readily recognized in 
this plot. For a mode at k = kj(a = 1) = 2.0/iMpc _1 , 
for example, it would continue to oscillate until now. Its 
oscillation 'frequency' is cu r ~ 0.76. Thus 2uj r ?] ~ 3.3 = 
0.52(27r), i.e., it would have oscillated for roughly half a 
cycle. By then its amplitude would be Landau-damped 
by more than one e-fold. In the last regime, the suppres- 
sion reaches an asymptotic value of ~ 70%. FigJ5] com- 
pares the conventional and the plasma-suppressed matter 
power spectra. We have replaced the shape parameter 
r = Sl m h in the BBKS formula fl2| for the transfer func- 
tion by the Eisenstein-Hu[l3| rescaled factor. Thus the 
baryon effects prior to the reionization epoch has been 
included in the transfer function, which is common to 
both cases. 
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FIG. 2: Matter power spectrum with and without the plasma 
effect. The conventions are the same as that in FigfT] 

By including the long-range residual electrostatic po- 



tential, we have shown that the growth of the baryon 
density contrast during the reionization epoch should be 
suppressed by the plasma collective effect in the form 
of ion-acoustic oscillations. Though subject to Landau 
damping, these oscillations happen to have both the pe- 
riod and the damping time comparable to the Hubble 
time. Thus there should exist certain imprints of these 
oscillations in the large scale structure formation history. 
In the linear regime of density contrast growth and the 
asymptotic limit where the oscillations damp away, the 
amount of suppression of the matter power spectrum due 
to the plasma effect reaches ~ 30%, which is sizable. In 
principle the plasma suppression is also effective in the 
nonlinear regime. It would be interesting to expand our 
treatment to the nonlinear regime of growth at late times 
(a — ► 1) and track its impact. 

One major effort in modern cosmology is to under- 
stand the cosmic composition through the determination 
of a set of cosmic parameters. In addition to the CMB 
fluctuations, large scale structure formation can help fur- 
ther constrain these parameters. In this regard, our pre- 
dicted plasma suppression of the power spectrum in the 
Lyman-a forest region may be relevant. We note that the 
observed power spectrum in the inter-galactic medium 
(IGM) region indeed appears to be systematically lower 
than that predicted by the conventional theory [lj|. It 
would be very interesting to confirm this. Aside from 
this issue, the evolution of the plasma-suppressed power 
spectrum should provide a unique window to reveal the 
detail history and dynamics of reionization. 
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